Wave packet autocorrelation functions for quantum hard-disk and hard-sphere 
billiards in the high-energy, diffraction regime 
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We consider the time evolution of a wave packet representing a quantum particle moving in a 
geometrically open billiard that consists of a number of fixed hard-disk or hard-sphere scatterers. 
Using the technique of multiple collision expansions we provide a first-principle analytical calculation 
of the time-dependent autocorrelation function for the wave packet in the high-energy diffraction 
regime, in which the particle's de Broglie wave length, while being small compared to the size 
of the scatterers, is large enough to prevent the formation of geometric shadow over distances 
of the order of the particle's free flight path. The hard-disk or hard-sphere scattering system 
must be sufficiently dilute in order for this high-energy diffraction regime to be achievable. Apart 
from the overall exponential decay, the autocorrelation function exhibits a generally complicated 
sequence of relatively strong peaks corresponding to partial revivals of the wave packet. Both the 
exponential decay (or escape) rate and the revival peak structure are predominantly determined by 
the underlying classical dynamics. A relation between the escape rate, and the Lyapunov exponents 
and Kolmogorov-Sinai entropy of the counterpart classical system, previously known for hard-disk 
billiards, is strengthened by generalization to three spatial dimensions. The results of the quantum 
mechanical calculation of the time-dependent autocorrelation function agree with predictions of the 
semiclassical periodic orbit theory. 



I. INTRODUCTION 

One of the main objectives of the field of quantum 
chaos is to study the quantum mechanics of systems that 
are chaotic in the classical limit, and to find what proper- 
ties of the classical systems appear crucial for the quan- 
tum description. There are some common approaches 
to the problem. One can explore the energy spectra of 
quantum analogs to classically chaotic, bounded systems 
H Q- LJ- lj- LJ > or investigate scattering resonances in the 
complex energy plane for unbounded quantum systems 
with chaotic classical repellers H, H, E3, El D3 • Inter- 
esting results relating the energy spectra of closed sys- 
tems to the scattering resonances of the complementary 
open systems (defined as the region "outside" the closed 
system) are also known 0]. Another approach focuses 
on the time-dependent description, which involves such 
quantities as t he q uantum state autocorrelation functions 
Q^I^ ITq |IT7llT^IT9l | and the Loschmidt echo, or fidelity 
[20l l2lL [23] ■ In this paper we continue the second ap- 
proach and address the time evolution of the autocorre- 
lation function for quantum particles traveling in hard- 
disk and hard-sphere billiards in two and three spatial 
dimensions respectively. 

We have previously studied the short time dynamics of 
a small Gaussian wave packet in arrays of hard-disk scat- 
terers 0]. The wave packet was considered to be much 
smaller than the disk radius, and its evolution was lim- 
ited to times shorter than the Ehrenfest time 0] , at which 
time the wave packet size becomes comparable with the 
scatterer size. The goal of the present paper is to extend 
the previous results by investigating the wave packet dy- 
namics in the limit of long times, much longer than the 



Ehrenfest time. Here, we show that substantial progress 
can be made for hard-ball, (hard disk or sphere) scatter- 
ing systems containing only a small number of scatterers, 
and for wave packets in the diffraction regime, to be fur- 
ther specified below. 

In a series of seminal papers 0, S 13 Gaspard and 
Rice studied classical, semiclassical and quantum prop- 
erties of a three-disk scattering system, in which a point 
particle moves in two dimensions in the presence of three 
fixed, impenetrable disk scatterers. The disks have the 
radius a, and their centers are placed at the vertices of 
an equilateral triangle of side R. Below we will refer to 
this scattering system as to the three-disk equilateral bil- 
liard to distinguish it from other three-disk billiards, in 
which the scatterer centers form isosceles or generic tri- 
angles. The latter scattering systems will be referred as 
to the isosceles billiard and the generic triangle billiard, 
respectively. The classical version of the three-disk sys- 
tem is known to have a fractal repeller with the positive 
Kolmogorov-Sinai entropy @ . The quantum dynamics of 
a particle in the three-disk billiard is largely determined 
by scattering resonances, whose distribution in the com- 
plex energy plane reflects the chaotic properties of the 
corresponding classical system H| ■ 

In this paper we address the time evolution of initially 
localized wave packets in such hard-disk scattering sys- 
tems as the two-disk billiard and the three-disk equilat- 
eral, isosceles and generic triangular billiards. We also 
extend our calculations to include other hard-sphere scat- 
tering systems in three spatial dimensions. As a tool for 
quantifying the wave packet dynamics we consider the 
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autocorrelation function C(t) defined as the overlap 

c(t)H(0 o |<MI 2 , (i) 

where |(/>o) is the initial quantum state of a billiard, and 
\4>t) is the quantum state obtained from the initial one by 
the time evolution through time t under the Hamiltonian 
H: 

\4>t)=e- im ' h \^). (2) 

We construct the quantum propagator, e ~ tHt / h ^ fo r a 
dilute hard-disk billiard using the multiple collision ex- 
pansion technique |U |24|. Then, we apply this prop- 
agator to an initially localized wave packet |</>o)j corre- 
sponding to a highly-energetic particle traveling through 
the system, and calculate the autocorrelation function 
C(t). The following two assumptions allow analytical 
treatment of the problem: (i) the scattering system is 
considered to be dilute, and (ii) the quantum particle is 
restricted to have a short de Broglie wave length. (We 
will give further details below.) These assumptions limit 
the validity of our results for the autocorrelation function 
to a certain time range. The latter can be quite long: it is 
bounded from below by the Ehrenfest time of the system, 
and typically extends over a large number of Ehrenfest 
times. Therefore, this work should be considered as an 
attempt to analytically describe the dynamics of initially 
localized wave packets over long times during which the 
quantum particle explores all the scattering system, and 
its wave function becomes spatially extended. 

It was pointed out by Gaspard et al. |9j that the time 
evolution of a wave packet in a hard-disk scattering sys- 
tem exhibits an overall exponential damping controlled 
by the location of the scattering resonances in the com- 
plex energy plane. As a result of this damping one ex- 
pects the wave packet autocorrelation function, C(i), to 
decay in an exponential-like manner as well. We ex- 
plicitly calculate the autocorrelation function for wave 
packets in various hard-disk billiards and investigate the 
structure of its time decay. We find that the decay con- 
sists of a sequence of sharp peaks with the exponentially 
decreasing envelope. The peaks reflect the phenomenon 
of the wave packet partial reconstruction due to interfer- 
ence of different classical periodic orbits in the billiard. 
The maxima of the peaks occur at times at which the 
classical particle, having the momentum equal to the av- 
erage momentum of the wave packet, would return to its 
initial phase-space point. Although the long-time auto- 
correlation function peaks have the same origin as the 
ones in the short-time autocorrelation function |16l Il9| , 
the peaks at longer times result from the constructive 
interference of a large number of classical trajectories 
whereas only a single classical path is responsible for the 
wave packet reconstruction during times shorter than the 
Ehrenfest time. 

The envelope of the autocorrelation function for a two- 
dimensional hard-disk billiard decays exponentially with 
time, e -7 *, with the rate 7 approximately equal to the 



difference of the mean positive Lyapunov exponent, A, 
and twice the Kolmogorov-Sinai (KS) entropy per unit 
time, /iKS, of the repeller of the corresponding classical 
system, i.e. 7 rs A — 2/ikS' This decay rate is nothing but 
the quantum escape rate originally obtained by Gaspard 
and Rice 0. Our generalization to three spatial dimen- 
sions shows that the envelope of the wave packet autocor- 
relation function also decays exponentially, e _73D *, with 
the decay rate 73D = Ai + A2 — 2haS) where the place 
of A in the expression for 7 is taken by the sum of the 
two mean positive Lyapunov exponents Ai and A2 of the 
classical hard-sphere system. 

The paper is organized as follows. In Section II we cal- 
culate the autocorrelation function, C(t), for wave pack- 
ets in two-disk and various three-disk scattering systems 
using the method of the multiple collision expansion. We 
show that the structure of C(t) is generally quite intri- 
cate, and we provide explanations of its main features. 
Section III presents an alternative description of the au- 
tocorrelation function decay based on the semiclassical 
arguments. While the semi-classical method described 
here does not provide as much detail about the autocor- 
relation functions as do the multiple scattering methods, 
it does provide information about the relative strengths 
of the peaks and their location. Moreover this method 
allows us to obtain information about more complicated 
systems that are difficult to treat using more exact meth- 
ods. For example, in this section we use the semi-classical 
method to treat generic three disk scatterers, and to 
generalize our methods to three spatial dimensions to 
describe the wave packet dynamics in two-, three- and 
four-sphere scattering systems. Our conclusions and dis- 
cussions are contained in the Chapter IV of the paper. 



II. HARD-DISK SCATTERING SYSTEMS 

A. General formulation 

We consider a particle of mass m placed among a col- 
lection of N fixed hard-disk scatterers of radius a cen- 
tered at position vectors Rj, with j = 1, 2, . . . , N. The 
Hamiltonian of the system can be written as 

N 

H = H + J2Vj, (3) 

where Hq is the free particle Hamiltonian, and the hard- 
disk scatterer potentials are given by 

V(r\-l +°° f ° r ' r ~ R J I < °' (A\ 
V ^ Y > - \ for jr-Rjl > a. W 

The time-domain propagator, G(t) = e - %Ht / h ^ satisfy- 
ing the Schrodinger equation with the Hamiltonian H, 
evolves an initial quantum state \<fio) in time in accor- 
dance with Eq. J2J. The corresponding energy-domain 
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propagator is obtained by means of the positive time 
Fourier transform 



G(E) 



1 

ih 



1 







E-H + iO 



The inverse transform is given by 

r+oo 



G(t) 



i 

2^ 



dEe- lEt ' h G{E), 



(5) 



(6) 



where we consider time t to be strictly positive. Although 
it is a formidable problem to construct the time-domain 
propagator directly, one can employ the method of mul- 
tiple collisions [23, [24J to obtain a series expansion for 
the energy-dependent propagator. 

The multiple collision expansion of the energy- 
dependent propagator represents G(E) as an infinite sum 
over possible collision sequences that the quantum parti- 
cle can undergo: 



G —Go + GqTjGo + GoTjGo?fcGo 

j j k=£j 

+ ^2 E E G TjG T k G TiG + . . . , 

j kjtj ¥fe 



(7) 



where Go = (E — H n + i0) _1 is the free particle prop- 
agator in the energy domain, and the binary collision 
operator, Tj, also known as the T-matrix, is defined by 
the multiple collision expansion, Eq. (J7J, written for a 
system containing only the j th scatterer, i.e. 



G, = — 

1 E 



1 



H - Vj + iO 



Go + GoTjGo. 



(8) 



Hereafter all G and T operators are given in the energy 
domain unless the time dependence is specified explic- 
itly. Here, Gj is the propagator for a particle moving in 
two dimensional space with only one scatterer located at 
point Rj. 

The momentum-space matrix elements of the binary 
collision operator for the hard-disk potential given by 
Eq. Q were calculated by Correia [2fj. For a particle 
of energy E = fi 2 K 2 /2m, where k is the magnitude of 
the particle's wave vector, the matrix element (k|Tj|k') 
of the binary collision operator, relating two generally 
different momentum states k and k', is given by 



(k|T J .|k')=2 7 ra^-e-( k - k ')^ £ 



Jl(8 k -8 k >) 



3l "' " 2m' 

1 — — OG 

7jf 2 [k'Jiik^Ji-xik'a) - *Ji_i(fca)Ji(fc' a)] 



k 2 - {k'Y 
+k'Ji{ka)Ji-i{k'a) - nJi{ka)Ji{k'a) 



(9) 



where (k,9k) and (k',8k>) are the polar coordinates of 
the wave vectors k and k' respectively, Ji is the Bessel 



function of the first kind of order I, hP^ is the Hankel 
function of the first kind of order I. Hereafter, the fol- 
lowing normalization conditions are adopted: (r|r') = 
5(r-r'), (k|k') = {2n) 2 S(k - k') and (r|k) = e jkr , so 
that the completeness relations read J dr|r)(r| = 1 and 
dk 

J |k)(k| = 1. Using Eq. © together with the ex- 

(■ ) 

pression for the free particle propagator matrix element, 



(klGolk') = 



2m 6{k-k r ) 



(10) 



Correia [25| has calculated the matrix element describing 
a sequence of n successive collisions of the particle with 
scatterers s,r,q, . . . ,p,j and i: 

(klTiGoTjGoTpGo . . . T q G T r G T s \k') 



= (-1)" ii e -»kR»+tk'R 3 

2m 



E 



Jlr\ Ka ) c ilr(6 ar -0rs) 



H, 



(i) 



(KRr S ) 



Ji s (k'a) ji s (e rs -e h ,-n/2) 
fli W (/so) 



(11) 



Here (Rij , 9ij ) are the polar coordinates of the scat- 
terer separation vectors Ry = R; — Rj, where i,j = 
1,2,3,. ..,N. 

Equation 10 together with l(TT|) allows one to calcu- 
late the energy-domain autocorrelation amplitude for an 
initial quantum state \(f>o), 



Sl(E) = (<p \G(E)\<Po) 

dk r dk' 



(2tt) 2 J (2tt) 



jxVO |k)<k|G(S)|k')(kVo) 



(12) 



The time-domain autocorrelation function, Eq. QJ, is 
then obtained from the energy-dependent overlap with 
the help of the inverse Fourier transform given by Eq. 10, 
namely 

c(t) = m)\\ (is) 

with the time-domain autocorrelation amplitude 

50 dE 



n(t) = {MG{t)\M = 



2tt 



e - tEt / R n{E). (14) 
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B. Wave packet 



tion over the k-space we obtain 



In order to calculate the autocorrelation function due 
to the propagator given by Eqs. Q and ijTT|l we consider 
a circular wave packet defined by 



Mr) 



e(a-|r-Ro|). 



(15) 



The wave packet represents a particle, with the average 
momentum Tiko, located at the position Ro. Here 6 is 
the Heaviside step function, and the real quantity a has 
an apparent meaning of the wave packet dispersion. Our 
choice of the wave packet is motivated by the demand to 
facilitate complex energy plane integration in the trans- 
formation of the autocorrelation overlap from energy to 
time domain, Eq. (|14|) . As we will see in Section III the 
time decay of the wave packet autocorrelation function 
does not depend significantly on the initial wave packet, 
as long as the latter is sufficiently localized in the coor- 
dinate and momentum space. 

Let us now fix the system of coordinates by imposing 
Ro = 0, so that the particle is initially located at the 
origin. The circular wave packet does not overlap with 
any of the disk scatterers if the following N conditions 
are satisfied: 



o\GoTiGoTjGoT p Go . . . T q GoT r GoT s Go\4>o) 
1 2m 
AiT 2 



,. . . ,l s — — oo 



JiA Ka ) 



jh(e 0t -e tj ) 



H^(na) 

Ji s i K a) e a e (0 re -e eO ) 



h[%,{kR s0 ) 



Xlr (K,k )e a '^- e "o-^) 



(19) 



The polar coordinates of the disk separation vectors 
R i = Ro-Ri and R s o = R s - Ro are given by (Roi, 9oi) 
and (RsO,Qso) respectively; the asterix denotes the com- 
plex conjugate. 



C. Diffraction regime approximation 



Rj>a + a, for j = 1, 2, 3, . . . , N. (16) 



The momentum representation of the wave packet is 
given by 



6 (k,k ) = 2^ 



Ji(|k-koH 



|k-k | 

+oo (17) 



with 



/cJ; + i(fccr)J ; (fcocr) - k J ; (ka) J/ +i (k a) 



k 2 — k 2 



(18) 

Here (ko, 9k ) are the polar coordinates of the wave vector 
ko. Appendix A provides the derivation of the above 
expansion. 

We are now in a position to calculate the part of the 
wave packet autocorrelation overlap due to a sequence of 
n successive collisions of the particle initially at Ro with 
scatterers s,r,q,... ,p,j and i. Performing the integra- 



The expression for the overlap in Eq. (|19J) is exact. We 
will now derive an approximation of this overlap for the 
case of a dilute scattering system. In the dilute hard-disk 
billiard, Rij 3> a, so that the high argument approxi- 
mation of the Hankel functions by exponentials can be 
used to greatly simplify the analysis of the autocorrela- 
tion function. 

Let us start with determining the angular momen- 
tum states dominating the autocorrelation function for 
a given value of k. The ratios inside the brackets in 
Eq. ((TT)|> . Ji{x)/h[ 1 \x), are small for I > x, since Ji{x) 
has its first maximum at x ~ I. (Figure ^ illustrates the 
latter argument for cases of / = 10 and I = 50.) Con- 
sequently, the main contribution to the multiple sum in 
Eq. I|19fl comes from terms with . . . ,l s running from 
— [na] to +[fta], and /,/' running from — [Ka] to +[kct], 
where the square brackets denote the integer part. At 
the same time, the large argument approximation of the 
Hankel function [HEll, 



holds for 



■ exp 



nl 
~2 



(20) 



(21) 



5 




x/l 



FIG. 1: The absolute value of the ratio Ji(x)/H^ (x) as a 
function of x/l for two cases: I = 10 and I — 50. The ratio 
significantly differs from zero only for x > I. 



Therefore, if uRij 3> (2na) 2 /2, and kRqi, kR s o 3> ( Ka + 
no) 2 /2 we can use this approximation in Eq. (|19|l to get 

1/2 



t'olGoTiGoTjGo . . . T r GoT s Go\4>o 



x (j>*(K,k ;6 0i - 9 ko )—=f K {9 a% - (%)- 



x f K 



i 

2-kk 




where 



MO) 



2 

7TIK 



1/2 +oo 

E 



(22) 



(23) 



is the scattering amplitude [24J describing scattering of 
a quantum particle of energy E = H 2 k 2 /2m from a hard 
disk of radius a at an angle 0. 

The approximation given by Eq. I|22|) is only valid for 
energies satisfying the conditions (see Eq. (|21|) and the 
discussion right below it) 



Ri 



2R 



0; 



2a 2 ' (a + af 



for i,j = l,2,3,...,JV. (24) 



These conditions bear a simple physical meaning. Sup- 
pose Rij, Roi ~ i? and a ~ a, then inequalities (|24|) can 
be written as i? 3> 2a/ a, with a = l/na. The latter has 
a meaning of the angle of diffraction of a wave with the 
wave length 1/k on an obstacle of size a, so that 2a/ a 
represents the estimate of the shadow depth, which is the 
largest distance over which the geometrical shadow can 
exist. Then, the conditions l|24|) simply mean that the 



scattering system is so dilute that the average separa- 
tion between scatterers is much greater than the shadow 
depth for the given particle's energy. This is equivalent 
to stating that Eq. 1221) is only valid in the diffraction 
regime, i.e. the diffraction effects prevail over the geomet- 
rical shadow, so that no disk scatterer can be screened 
from the particle by other disks. In the case of a suf- 
ficiently dilute scattering system, R 3> a, the inequali- 
ties (|24() are satisfied for a significant range of energies, 
so that Eq. ijl'l'l) provides a good approximation of the au- 
tocorrelation function for wave packets confined to this 
energy range. 

Equation (1221 has an apparent structure. The initial 
wave packet \4>q) is propagated through a sequence of free 
flight and collision events, each of which contributes by 
a corresponding product term to the expression for the 
autocorrelation overlap. Indeed, the latter is multiplied 
by e lKi? /Vi? each time the particle experiences a free 
flight through some length R, and by f K (8) each time 
the particle scatters off a disk at an angle 9. The same 
wave function construction algorithm was earlier used in 
reference for semiclassical quantization or the three- 
disk scattering problem. Nevertheless, it is important 
to note that, at least for the autocorrelation function 
calculation, this method fails in the true semiclassical 
limit, k — > oo, due to violation of the conditions given 
by 124|) . and is legitimate only in the diffraction regime 
approximation defined by 124(1 . 

A closer examination of Eq. 1221 shows that only such 
initial wave packets that represent particles moving in 
the vicinity of system's classical periodic orbits can ex- 
hibit a significantly non-zero autocorrelation function for 
times corresponding to a number of collision events of 
the counterpart classical particle. To see that, suppose 
that the wave packet \<f>o) given by Eq. (|17|l is well lo- 
calized in momentum space, i.e. the particle's dc Broglic 
wave length X = 1/fco <C o . Then, the overlap given 
by Eq. ll2*2*|l is negligible unless 9 0i ss 9 s0 m 9k . Indeed, 
4>o(k, ho; 6), if considered as a function of the angle 9 be- 
tween the wave vectors k and ko, is sharply peaked at 
9 = 0, and rapidly vanishes as k turns away from k . 
This means that the autocorrelation overlap (<f>o\G\<f)o) 
gets a significant scattering contribution in addition to 
its free-streaming part, (<f>o\Go\(f>o) , (which is negligible 
for dilute systems) only if the wave packet is initially lo- 
cated on and moves along a line connecting centers of any 
two disks in the scattering system. This happens because 
the reflection wave produced by a disk at the last scat- 
tering interferes destructively with the initial wave unless 
the two waves have their wave vectors pointing almost in 
the same direction. Hence, one can expect substantially 
non-zero values of the autocorrelation function only for 
such initial wave packets that represent classical parti- 
cles moving along lines connecting scatterer centers, and 
therefore traveling in the vicinity of unstable periodic or- 
bits of the hard-disk scattering system. 

We will now use the propagator given by Eqs. and 
(12211 to calculate the autocorrelation function for the cir- 
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cular wave packet defined in Eq. (|17|) in three different 
scattering systems: two-disk, three-disk equilateral and 
three-disk isosceles billiards. The autocorrelation overlap 
in the energy domain can be written as 



(25) 



where £lo(E) = (0o|Go|0o) is the free flight part of the 
overlap, and the scattering part, Cls(E), is determined 
by all possible collision events: 

3 

+ Y,T,^o\GoT 3 G T k G \cbo) 

3 kjtj 

3 k=£j lyik 

(26) 

The series in Eq. 1|26[1 can be summed explicitly for the 
three above-mentioned billiard systems using the diffrac- 
tion regime approximation, Eq. I|24fl . together with the 
assumption of that the initital wave packet has suffi- 
ciently high energy, as will be described below. 



D. Two-disk billiard 



where the diffraction regime approximation, Eq. (1221) , has 
been used. Here, r\ and r 2 are the distances separating 
the center of the wave packet and the centers of disks "1" 
and "2" respectively; n + r 2 = R. Substituting Eq. J27J| 
into Eq. (|26|l . while neglecting all other scattering se- 
quences, one obtains a geometric series that sums to 



n s (E) 



iR 



Ti \2nKrir 2 



1/2 



|<M K > k ; 



1 



„ikR 



(28) 



Equation (|24|l along with the assumption r\ ~ r 2 ~ R/2 
shows that the validity of the result predicted by Eq. I|28|) 
is limited to the energies satisfying n <C R/2a 2 , R/(a + 
a) 2 . Now, if the wave packet is sufficiently localized in 
the momentum space, i.e. the dc Broglie wave length 



A = l/fco < <r, 



(29) 



and if k < R/2a 2 , R/(a + <r) 2 , then Eq. can be used 
to calculate the time domain autocorrelation function for 
the wave packet. 

Finally we make an assumption that the quantum par- 
ticle is highly energetic, so that its de Broglie wave length 
is much smaller than the scatterer size, i.e. A <C a. 
This approximation allows us to use the semiclassical 
(WKB) expression for the hard-disk scattering ampli- 
tude, namely 




FIG. 2: Two-disk billiard. The circular wave packet is initially 
located on the classical periodic orbit distance ri away from 
the center of disk "1" , and distance r 2 away from the center 
of disk "2" , with ri + r 2 = R. 



The simplest of all hard disk systems is the two-disk 
billiard, which consists of two hard disks, "I" and "2", 
of radius a with the center-to-center separation R, see 
fig. [3 Following the discussion above we place the initial 
wave packet on the line connecting the disk centers with 
its wave vector ko pointing along this line in order for 
the sum in Eq. (|26|l to have significant, non- vanishing 
terms. For the initial condition shown in fig. |3 these 
non- vanishing terms are 

(4>o\GoTiGoT 2 Go ■ ■ ■ TiG T 2 Go\4 l o) 

g£rer*i giizR ^ikR 

X -— f K (Tc)—— .f K ( K )—— 

\AT Jr Jr 



' /kW 7I /kW ^ 



He) 



|sin(0/2)| e -2i«a|Bin(e/2)|_ ( 30 ) 



Equation (|30l) is known to be a good approximation of 
the exact scattering amplitude for sufficiently large scat- 
tering angles |28j |. Then, combining the diffraction ap- 
proximation, Eq. (|24|l . with the short de Broglie wave- 
length approximation, A <C a, we arrive at the following 
condition: 



2a 2 (a + cr) 2 
~R 



R 



« A < a. 



(31) 



We refer to this inequality as to the condition of the high- 
energy diffraction regime. In the rest of this paper we 
assume that both the momentum space localization con- 
dition, Eq. I|29l) . and the high-energy diffraction regime 
condition, Eq. 1|) . are satisfied. 

The calculation of the autocorrelation function in the 
time domain, C(t), is carried out in accordance with 
Eqs. (H3Jl, G2}, (EH and (SgJ. The scattering part of 
the time-domain autocorrelation amplitude reads 



ils(t) = / annexp 

2tt m Jy 



(27) 



where ils(E) is given by Eq. 1)28(1 ■ and contour T in the 
complex K-plane follows the imaginary axis from —ioo to 
0, then turns at the right angle, and proceeds to +oo 
along the real axis. 
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Careful analytical calculation of a range of times, for 
which Eq. with Q S (E) given by Eq. i|28|) . yields 

accurate predictions, is a formidable problem. Neverthe- 
less, we can use simple physical arguments to roughly 
estimate this time range. The method that we used to 
calculate the energy- dependent scattering part of the au- 
tocorrelation function, Clg(E), relies on the assumption 
of the wave packet diffraction. The wave packet must ex- 
plore the scattering system for the diffraction effects to 
take place. An estimate of the time needed for the parti- 
cle to reach the first scatterer is ~ r 2 /v ~ R/v, where 
T2 is the distance between the initial location of the par- 
ticle and the first scatterer it collides with, see fig. [3 and 
v = hko/m is the average velocity of the wave packet. 
Time t e also gives an estimate of the Ehrenfest time for 
the system: for times shorter than tE the wave packet 
evolution is dominated by the free particle Hamiltonian, 
and therefore is classical-like, while particle's propaga- 
tion is diffractive and substantially non-classical for times 
beyond tE- To estimate the upper bound of the applica- 
bility time-range we note that the WKB expression for 
the scattering amplitude, Eq. (|3U|) . breaks down at low 
energies E = h 2 K 2 /2m with k ~ I/a. The momentum 
corresponding to these energies is Hk ~ (X/a)hko, and 
the corresponding velocity is v K ~ (X/a)v. The contri- 
bution of these low energy modes of the particle to the 
autocorrelation function become significant after the long 
wavelength part of the wave packet explores the scatter- 
ing system, i.e. after times i max ~ R/v K ~ (a/X)R/v cor- 
responding to a I X particle-disk collisions. Since a ^S> X 
this number of collisions can be quite large. Summa- 
rizing, we find that the applicability time-range of our 
analysis of the time-dependent autocorrelation function 
is estimated by 



f < vt/R<a/X. 



(33) 



In order to evaluate the integral in Eq. ilML'l) one can 
show that for t > the contour can be closed along the 
infinite quarter-circle k = Kooe 17 , with — > +oo, and 
the angle 7 decreasing from to — 7r/2 |4l|, so that the 
value of the integral is determined by poles of Q,§{E) in 
the fourth quadrant of the complex k plane: 



nn i 2R 

Kn = 111 , 

R 2R a ' 



for n = 1, 2, 3, 



(34) 



The semiclassical approximation for the scattering ampli- 
tude, Eq. (|3U[) . was used to find zeros of the denominator 
of fig (£7), so that Eq. (|34|l correctly locates the poles 
right below the region on real k axis where 4>o( K , ko] 0) is 
localized. Then, calculating the residues corresponding 
to the poles, we get 



fi s (t) S3 - (2mR ri r 2 ) 



n + [R/a] 

x E 

n ~ no — [R/ or] 



-1/2 



H^\(j) {K ni k \Q)\ z exp ( -i^-n 2 n ) , 

2m 



Tit 



where ttuq/R = ko, and the square brackets in the sum- 
mation limits represent the integer part. Equation 135|) 
is expected to hold within the time range given by 

The free streaming part of the autocorrelation function 
overlap, fio(i) = (<j>o\Gn(t)\(j) ) , is calculated explicitly 
for the free particle propagator. It can be shown that its 
contribution to the full autocorrelation function C(t) is 
negligible for dilute billiard systems j^- Therefore, the 
autocorrelation function for long times is entirely deter- 
mined by the scattering events, so that C(t) » |fig(<)| 2 . 

The main features of the time-domain autocorrelation 
function C (t) can be deduced by considering only a small 
number of poles with n = no + h, where h is sufficiently 
small, so that the pre-exponential function in Eq. I|35|) 
stays approximately constant. The contribution due to 
these poles is 



n s(*) ~ E exp 

n 

-iBot/n 



fit 1 7; n 

-l- I Ko 



2m 



i 2R 
~R~2R n ~ 



exp (-hVtWe-* 1 * 



(36) 



where Eq = fi 2 k 2 /2m is the average energy of the wave 
packet, v = hko/m is its average velocity, and 



it a 



(37) 



is the classical Lyapunov exponent of the two-disk pe- 
riodic orbit in the limit R ^> a, e.g. see reference |30j . 
Equation (|36|l shows that (i) the envelope of the scat- 
tering part of the autocorrelation function decays expo- 
nentially with time, C(t) ~ e~ A< >4 , with the decay rate 
given by the classical Lyapunov exponent of the system, 
and (ii) strong interference peaks occur in C(t) at times 
t multiple to the period of the classical periodic orbit, i.e. 
when vtj R is an even integer. The peaks of the autocor- 
relation function correspond to partial reconstruction of 
the wave packet at times at which the counterpart classi- 
cal particle returns to its initial point in the phase space. 

Figure |3| shows the autocorrelation function, C(t), 
computed using Eq. with the wave packet given 

by Eq. (|I7|I . for the two-disk billiard with the following 
parameters: a = 1 and R — 10 4 . The circular wave 
packet of size a = 1 is initially placed in the middle be- 
tween the disks, r\ = r 2 = R/2, and its de Broglie wave 
length A = I0~ 2 . Conditions and IpHJl are satis- 
fied by this system. The summation in Eq. I|35(l includes 



20,001 poles. The solid line shows e 



-A< 2 >t 



decay. The fig- 



(35) 



ure shows the decay for times t greater than the Ehrenfest 
time tE ~ R/2v. 

As mentioned above, fig. |21 exhibits the wave packet 
reconstruction (revival) peaks together with the expo- 
nential decay of the autocorrelation function envelope. 
Another distinct feature of the decay is the broadening 
of the peaks in the course of time. These phenomena 
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vt/R 



FIG. 3: Wave packet autocorrelation function C(t) vs. vt/R 
for the two-disk billiard. Parameters of the system are as 
follows: a = a = 1, R = 10 4 , n = r 2 = i?/2 and A = 10~ 2 . 
The straight line shows exponential decay with the rate given 
by the classical two-disk Lyapunov exponent A' 2 '. The decay 
is shown for times t greater than the Ehrenfest time tE ~ 
R/2v. 

have a simple physical explanation. The direction along 
the two-disk periodic orbit is neutrally stable, i.e. a per- 
turbation of the initial phase-space point (of a classi- 
cal particle on the periodic orbit) along this direction 
growths at most linearly with time. On the other hand, 
the direction perpendicular to the periodic orbit is ex- 
ponentially unstable, with the Lyapunov exponent A' 2 ' 
playing the role of the classical instability rate. There- 
fore, the wave packet probability density dies out lin- 
early with time in the direction along the periodic orbit, 
whereas it decreases exponentially in the perpendicular 
direction (S^. The spreading of the wave packet in the 
direction along the periodic orbit (and therefore along 
its average velocity) yields prolongation of time intervals 
during which \<j> t ) significantly overlaps with the initial 
state |0o) ■ This prolongation amounts to broadening of 
the revival peaks of the autocorrelation function. On the 
other hand, the wave packet spreading in the unstable 
direction is responsible for the exponential decay of the 
autocorrelation function envelope. 

E. Three-disk systems 

We will now calculate the autocorrelation function de- 
cay for different scattering systems, namely three-disk 
billiards, in which a particle moves in two dimensional 
space with three fixed hard-disk scatterers. Here we re- 
strict ourselves only to such three-disk billiards for which 
the centers of the disks constitute vertices of an equilat- 
eral or of an isosceles triangle. The corresponding scat- 
tering systems are then referred to as the three-disk equi- 



lateral and isosceles billiards respectively. The discussion 
of the generic three-disk billiard, in which all three sides 
of the triangle are different, is left for the Section III. 

1. Three-disk equilateral billiard 




aR aR 



2 



'l '2 

FIG. 4: Three-disk isosceles billiard. The circular wave packet 
is initially located distance r\ away from disk "1", and dis- 
tance ri away from disk "2" , with r\ + r 2 = R. Disk "3" is 
distance aR away from disks "1" and "2". The "equilateral" 
billiard case corresponds to a = 1. 

We consider a geometrically open billiard consisting of 
three hard disks of radius a centered at the vertices of an 
isosceles triangle with one side of length R and two other 
sides of length aR, see fig. 0J Here we focus on the case 
a = 1, which allows complete analytical treatment. The 
following section is devoted to the case a ^ 1, for which 
a substantial understanding can be achieved in the limit 
a > 1. 

As in the case of the two-disk billiard we initially put 
the wave packet of size a on the line connecting the cen- 
ter of two disks, labeled by "1" and "2", with center-to- 
center separation R, see fig. 0] Distances r\ and r 2 sep- 
arating the wave packet center and the centers of disks 
"1" and "2" respectively satisfy the apparent condition: 
r*i + r2 = R. The average wave vector of the wave packet 
ko is pointing along the line connecting the centers of 
disks "1" and "2" as shown in fig. 

Our calculation of the wave packet autocorrelation 
function employs the transition matrix method first ap- 
plied by Gaspard and Rice yj for calculation of scattering 
resonances in three-disk scattering systems. Transition 
matrices, and related to them monodromy matrices, have 
also been used in analysis of different systems by Bogo- 
molny 31], and Agam and Fishman 32]. We fist spell 
out the multiple collision expansion given by Eq. i|2rj|l for 
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the case of a three-disk billiard: 



&s(E) = 51 (4>o\GoTiGoTiGoTj . . . G T 2 Ga\cpo/, 

n=2 path(rt) 

(38) 

where the second sum runs over all possible paths con- 
sisting of n binary collision events with the first collision 
taking place at disk "2" and the n th collision at disk "1" . 
Every term in this double sum is evaluated using the 
diffraction regime approximation, Eq. (|22|l . Following 
we construct a 6 x 6 matrix Q describing a transition 
in 6-dimensional space composed of directions (1 — *■ 2), 
(1 -> 3), (2 -► 1), (2 -> 3), (3 -> 1) and (3 -> 2), due to a 
single scattering event, 
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(39) 



where 



X = f K (ir)— and W = /„(2tt/3)-=. (40) 



i? 



i? 



Here 7r and 27r/3 are the turning angles for a classical 
particle bouncing among three disk of radius a placed in 
the vertices of an equilateral triangle with side R ^> a. 
The second sum in Eq. I|38|l is then given by the one-one 
element of the matrix Q": 



(^o\G TiG TiGoTj . . . G T 2 G \^o) 

path(n) 

- _™ ( lR X V2 

ft 2 \2irnr\r2 



(41) 



|0o(^fco;O)| 2 (Q™) 11 . 



Substituting Eq. 141|l into Eq. I|38|l , and taking advantage 
of the equality Y^,n=2 Q" = Q 2 (1 Q) > we obtain 



^(^) = 5f^rV /2 |^( Kl fco;0)| 2 



x l- 



1/6 1/6 



1-W + X 1-W-X 

(2 + W)/3 



(42) 



1 + W + W 2 - X 2 



where j 



1,2,3,4, and 
£i 



1 



6 



(l/2) 1/2 + (3V2/ 4 ) 1/2 ' 

/2 _ 3 3/2\V2 _ 3 i/ 4 



( 4 3/ 



1 



(3/4) 1/2 
1/2 



( 4 3/2_ 3 3/2)V^ + 3 l/4 



(44) 



£4 = 



1 - (3/4) 1/2 
1 



(1/2) 



1/2 



(3V2/4) 



1/2 



Equations l|43|l and (|44|l . which locate the scattering reso- 
nances of the three-disk system, were originally obtained 
by Gaspard and Rice 0. It is interesting to note that 
K n2 and k„3, being the simple poles of the autocorre- 
lation amplitude Cls(E), appear as double poles in the 
three-disk scattering matrix 0. The three-disk scatter- 
ing resonances have been also calculated by Cvitanovic 
and Eckhardt 0] b y m eans of the fundamental cycle 
expansion technique |33j | . 

The time-domain scattering part of the autocorrela- 
tion function, f2s(^) 5 is determined following the proce- 
dure used for the analysis of the two-disk billiard sys- 
tem. As before we only consider the poles lying under 
the region on the real /t-axis on which the wave packet is 
mainly concentrated, i.e. n E (n — [R/2a], no + [R/2a]), 
with 2nno/R = fco and the square brackets denoting the 
integer part. Calculation of residues of Sls(E) is straight- 
forward. The result is given by 



n s (t) w - (2mRnr 2 ) 
6 



-1/2 



n + [R/a] 
~ no-[R/2o\ 



/ «^I|'/ , o(K„i,fco;l 



(45) 



+ 2y^|0 o (^2,A:o;O)| e~ l ^ K ^ 

+ 2^i \(j> {K n3 , fc ; 0)| 2 e - l ^ K ^ 

+ \A^4 |0o(K„ 4 ,fco;O)| 2 e" 4 ^r K ^| . 

In order to predict the main features of the decay one 
can consider only poles in a small vicinity of the peak of 
the wave function, n = n + n with n -C n . This yields 



+ 2e -72t/2 + e -74t/2 



iE t/h 



(46) 



As in the two-disk billiard case the poles of Qs(E) lo- 
cated in the fourth quadrant of the complex k plane de- 
termine the time-domain autocorrelation function. Fol- 
lowing we define £ • 



-V^/R 



to find that in the 



limit R 3> a the poles of S (E) are given by 
27m + 7T + arg£j 



R 



2R a 



(43) 



where, as before, Eq is the average energy of the 
wave packet, v is its average velocity, and jj — 
(v/R)\n(R\^\ 2 /a) with j = 1,2,3,4. The slowest de- 
cay rate, 71, playing the role of the decay rate of the 
autocorrelation function envelope, can be written as 

4R v , OMR 

1"" 



51 =71 = ^1° — ":„ l2 «~ln-^--. (47) 



[2V2 + 31/4] 



w 1 
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Once again neglecting the free-streaming amplitude 
fioW we calculate the autocorrelation function as C(t) rs 
|ils(^)| 2 - As m the two-disk billiard case C(t) exhibits a 
sequence of strong wave packet revival peaks correspond- 
ing to phase space returns of the counterpart classical 
particle, see fig. El Taking into account that |£i| ~ 0.73, 
| £2 1 « 1.34, |f 3 | w 11.16 and |f 4 | ~ 20.38 one can notice 
that the peaks occur whenever vt/R is an integer greater 
than one, see fig. [SJ At t = R/v a part of the wave 
packet reflected by disk "2" overlaps with the initial wave 
packet, but this overlap leads to destructive interference 
since the wave vectors of the two waves have opposite 
directions. This shows up as the absence of the revival 
peak at vt/R = 1, and appears mathematically as partial 
cancellation of the expression within the square brackets 
in Eq. (|4(j|l . ft is evident from Eq. (|46[) together with the 
inequality |£i|,|£2| <C I&IjI&I that, as pointed out by 
Gaspard and Rice Q , the lines of poles corresponding to 
£3 and £4 are screened by the other two lines of poles and 
have no affect on the wave packet dynamics. One can also 
show that after only two collisions the RHS of Eq. I|4tj|) 
becomes totally dominated by the first exponential term 
within the square brackets resulting in the exponential 
decay of the autocorrelation function, C{t) ~ e -7<3> *. 



1, R = 10 4 , n 



'2 



R/2 




FIG. 5: The autocorrelation function vs. vt/R for the three- 
disk equilateral billiard. Parameters of the system are as fol- 
lows: a = a = 1, R = 10 4 , n = r 2 = R/2 and A = 10~ 2 . 
The straight line shows exponential decay with the rate given 
by 7 ( - 3 - > . The decay is shown for times t greater than the 
Ehrenfest time tE ~ R/2v. 

Figure [S] shows the decay of the time-dependent auto- 
correlation function for the three-disk equilateral billiard 
studied above. The parameters of the system are chosen 
to be identical with once used in the case of the two- 



disk billiard: a = c 

and A = 10~ 2 . The function C(t) for a time interval 
comprising 10 collision events was calculated by comput- 
ing the sum in Eq. H45|) over the total of 40,006 poles. 
The straight line shows e~ 7 4 decay, with 7^ given by 
Eq. (|47|l . It is interesting to note how small the magni- 
tude of the autocorrelation function becomes after only 
a few particle-disk collisions. After the time correspond- 
ing to ten bounces of the classical particle, the return 
probability drops down to a value below 10 -40 implying 
practical orthogonality of the initial and final states of 
the quantum particle. It is the consequence of the huge 
scale difference in the billiard considered here. The phe- 
nomenon of the wave packet partial reconstruction can 
be enhanced by decreasing the ratio R/a. 



2. Three-disk isosceles billiard 

We now address a more general three-disk scattering 
system, in which the scatterers are located at the vertices 
of an isosceles triangle, as shown in fig. 01 with a / 1. 
Once again the circular wave packet is initially placed 
between disks "1" and "2", which have center-to-center 
separation R, while disk "3" is distance aR away form 
them, see fig. 0] 

The single collision transition matrix Q, which in the 
case of an equilateral billiard was given by Eq. (|39|) , now 
reads 



1-2 1-3 2 1 2-3 3 1 3-2 
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(48) 



with 



X x = Un) 



R 



X 2 = / k (tt) 



Wx = f K (7T - fa) 



a itzR 



R 



w 2 = f K {-K-<t>x)- 

■IKO.R 



(49) 



W 3 = U(*-<hY 



Here, <fix and </> 3 are angles of the triangle corresponding 
to vertices "1" and "3" respectively; the angles are func- 
tions of a, and satisfy the obvious relation: 2(f>x + 4>3 = n. 

Following the technique used above, we calculate the 
one-one element of the matrix Q(l — Q) 1 to obtain the 
energy dependent autocorrelation function: 
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n s (E) 



rn 

~2 



iR 



h 2 \2irKrir 2 

xKi 



1/2 



|^o(k, fc o ;0)|' 

[2X 1 X 2 + W 3 - 



Xl(2X x X 2 - VKiVK 2 )] W X W 2 - {X ± X 2 



W{W 2 ) 2 W? 



(50) 



[1 - XI + (JCxXs - W 1 W 2 )W 3 \ Z - - *|) + XaOW 2 + W 3 



As above, we introduce £ = —Wa/Re 1 . The poles 
of Sls(-E), i.e. zeros of the denominator in Eq. (|50|l . are 
given by solutions of the following two polynomial equa- 
tions: 



2+ V2- 



where 




'■2a 



LU- 



-r) 



R 



2 + %/2£ 
(51) 



2a 



u± = 



1± 1- 



4a 2 



1/4- 



(52) 



In the limit of a — > +oo we have lo + — ► a and w_ — ► 2 5 a 3 . 

For the sake of clarity of the following analysis we as- 
sume 2a to be an integer. Thus, for R 3> a and large a 
one can find approximate solutions to Eqs. (|51(l : 
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FIG. 6: (a) First two "bands" of poles for the case of a = 5/2 
and R/a = 10 4 ; (b) magnification of the first band; (c) mag- 
nification of the second band. Dots correspond to (exact) val- 
ues of the poles numerically computed directly from Eq. 1511 , 
while crosses show the same poles approximated by Eq. (1531 . 
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1-1/q 
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(0) _ 1/2q 
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1-1/q 



c ,i7r(l+p/a) 




(53) 

where p — 1,2,..., 2a. Each of these 4a + 2 values of £ 
defines a line of poles of fig(£0 m the complex K-plane 
in accordance with 



27m + 7r + arg £ Pj 



(±) _ 27m + 7r + arg£ 

Si, 3 — 



(±) 



R 



2R 



2R a 



with j — 1, 2; 

(54) 



The poles given by Eqs. (|54|l come in three "bands" . 
The first and the second bands of poles, corresponding 
to £ p i and £ Pi 2 respectively, are essential for the wave 

packet dynamics, while the third band, given by £3 , is 
completely screened by the first two due to the relation 
|£i| < I £2 1 I £3 1- Figure HJ1 shows the first two bands 
of poles for the three-disk isosceles billiard with a = 5/2 
and R/a = 10 4 . Figure displays the bands over some 
interval of real K-axis, while fig. [SJd and fig. magnify 
the first and the second bands respectively. The dots in 
the last two figures represent poles numerically computed 
directly from Eqs. i|5"T)l. and thus should be thought as of 
"exact" poles, while crosses are the poles approximated 
by Eq. 1|53[1 . In most of the cases the dots and the crosses 
fall on top of each other, so that one can conclude that 
Eqs. (|53|l accurately locate poles of the autocorrelation 
function. 

As we saw earlier, the size of the gap separating the 
poles and the real K-axis determines time decay of the 
envelope of the autocorrelation function. Thus, the auto- 
correlation function envelope decay, C(t) ~ exp(— 7^ t), 
for the isosceles three-disk billiard is governed by the rate 
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In the limit a 3> 1 the decay rate becomes 7a 
(v/aR) \n(aR/a). 



(55) 

(3) _ 




FIG. 7: The autocorrelation function as a function of time for 
the isosceles three-disk billiard with a = 5/2. Disks "1" and 
"2" are separated by distance R — 10 4 a. The wave packet 
of the de Broglie wavelength A = 10~ 2 a is initially located as 
shown in fig.^Jwith ri = T2 = i?/2. The dotted and the solid 

straight lines represent e _A< '* and e _7 ° * decays respectively. 

The time-domain autocorrelation function, C(i) ~ 
|Sls(i)| 2 , is now calculated using Eq. (23, with fls(E) 
given by Eq (| 5 0|> . As before, the complex K-plane con- 
tour integration is performed by calculating the residues 
of Qs (E) and computing the sum over the poles given by 
Eq. (52} with the index n £ (no - [R/2a],n + [R/2a]). 
Here 2-ktiqIR = ko and the square brackets denotes the 
integer part. The autocorrelation function C(t) for the 
three-disk isosceles billiard is shown in fig. [7| for the case 
of a = 5/2. The system is characterized by R/a = 10 4 ; 
the wave packet of the de Broglie wavelength A = 10~ 2 a 
is initially placed in the middle between disks "1" and 
"2", see fig. H with r 1 = r 2 = R/2. Figure □ confirms 
our earlier predictions that the overall envelope of the 
autocorrelation function decays as e _7 ° * with the decay 
exponent given by Eq. (|55|l . The trend of this expo- 
nential decay is presented by the solid lines in the figure. 



The dotted trend line corresponds to e 



decay, with 



the rate A^ 2 \ Eq. (|3*7|) . being the two-disk Lyapunov ex- 
ponent of the unstable periodic orbit trapped between 
disks "1" and "2". 

At first glance the autocorrelation function C(t) ap- 
pears merely as some complicated sequence of decaying 



peaks. Nevertheless simple physics underlies its struc- 
ture. The peaks of the autocorrelation function, or the 
wave packet partial revivals, occur at instants of time at 
which the counterpart classical particle returns to its ini- 
tial point in phase space. Then, relatively larger peaks 
result from the phase space trajectories with the smaller 
number of collisions per unit time interval. It is because 
at every collision event a dominant part of particle's prob- 
ability density completely escapes the billiard, and only 
a tiny part of this density proceeds to the next collision 
in order to eventually contribute to the autocorrelation 
function. In other words, phase space periodic orbits 
with longer mean free paths result in stronger reconstruc- 
tion peaks. In the isosceles three-disk billiard with a > 1 
the long free flight path trajectories are the ones that pass 
through disk "3" every second collision, see fig. 0] and re- 
turn to the initial point at times t n = (2cm + l)R/v with 
n = 1,2,.... Thus, C(t) exhibits relatively strong peaks 
at times t n corresponding to the scattering sequences 
"132", "13232", "13132" etc. During time intervals be- 
tween the large peaks, i.e. i n _i < t < t n , smaller wave 
packet reconstruction peaks occur due to periodic orbits 
with shorter mean free paths, which are determined by 
trajectories bouncing mostly between disks "1" and "2". 
This is how the two-disk Lyapunov exponent A^ 2 ) enters 
the description of C(t) for the three-disk billiard. Thus, 
for times t < t\ only the two-disk collision sequences 
"12" , "1212" etc., contribute to the autocorrelation func- 
tion, resulting in the e~ A< * decay. 

Another distinctive feature of fig. [7| is the absence of 
revival peaks for times R/v, 3R/v and 5R/v. That 
is because for a classical particle moving with velocity 
v = hko/m there are no phase space periodic orbits cor- 
responding to these times. We have already encountered 
same phenomenon while discussing the autocorrelation 
function decay in the three-disk equilateral billiard. 

Now it is also apparent how the two-disk decay e _A 2 * 
is recovered if disk "3" is removed to infinity, see fig. 01 
In the limit a — > oo the time t\ of the first large re- 
construction peak goes to infinity, resulting in the two- 
disk billiard decay of the autocorrelation function for 
t < t\ — ► oo. 

Before closing this section we give a brief summary of 
the results obtained. We have used the multiple collision 
expansion technique to provide a detailed first-principle 
calculation of the time dependent autocorrelation func- 
tion, C(i), for quantum wave packets moving in billiards 
composed of two and three hard disk scatterers. By an- 
alytically constructing C(t) for three-disk isosceles bil- 
liards we have broadened the class of three-disk systems 
earlier treated by similar methods 0- The applicability 
limits of these methods were shown to be given by the 
condition of the high-energy diffraction regime, Eq. 131fl . 
We found the decay of the autocorrelation function for 
the case of a three-disk equilateral billiard to be mainly 
exponential apart from a regular sequence of wave packet 
reconstruction peaks of equal relative strength. On the 
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other hand, C{t) decays non-uniformly in the case of a 
three-disk isosceles system, and features different decay 
rates at different time scales. Thus, at times shorter that 
(2a + l)R/v, see fig. 0] the decay is entirely determined 
by the Lyapunov exponent of the shortest two-disk clas- 
sical periodic orbit, while at long times, the overall en- 
velope decays at a slower rate given by Eq. I|55|l . A well 
pronounced structure of revival peaks is again observed 
in the autocorrelation function of the three-disk isosceles 
billiard. As we will show in the Section III the relative 
strengths of the peaks are determined by the number 
of interfering periodic orbits of the counterpart classical 
system. 



III. SEMICLASSICAL DESCRIPTION 

In this section we present a simple method for pre- 
dicting relative strength of the peaks of the wave packet 
autocorrelation function C(t), see Eq. based on the 
semiclassical Van Vleck propagator [34|. The semiclas- 
sical analysis of the autocorrelation function was earlier 
performed by Tomsovic and Heller for the case of the 
stadium billiard (35[. We start with applying the semi- 
classical method for two- and three-disk billiards studied 
in the previous section, and then, use it to investigate 
such more complicated scattering systems as the three- 
disk generic billiard and the two-, three- and four-sphere 
billiards in three spatial dimensions. 



3. Semiclassical approach 

In the limit of short de Boglie wavelengths the time 
evolution of a quantum state in a hard-disk or hard- 
sphere billiard can be described by the semiclassical Van 
Vleck propagator 



C? sc (r,r / ;t) = (r|G sc (t)|r / ) 



2irih 



d/2 



2^,DJ exp I i- 

7 



(56) 

The summation in this expression goes over all classi- 
cal paths 7 connecting points r' and r (in ei-dimensional 
space) in time t. S 7 (r,r';i) represents the classical ac- 
tion along the path 7, and is an index equal to twice 
the number of collisions of the particle with hard scat- 
terers during time t 0. The Van Vleck determinant, 
D 7 = I det(— d 2 Sj/drdr')\, corresponds, up to an appro- 
priate normalization factor, to the classical probability of 
the path 7 [36|. 

The autocorrelation overlap due to the propagator 
given by Eq. (|56|l and the an initial quantum state |</f>o) 
can be written as 



o\G sc (t)\<t>o) = / dr / dr'^(r)G sc (r,r';O0o(r'), 



where asterix denotes complex conjugate. Let us now 
assume that the wave function <fio( r ) is localized about 
point Rq and has an average momentum ?iko, 



ik r 



= \Mr)\e* k0 



(58) 



Since both wave functions in the integrand on the right 
hand side of Eq. (|57|l are localized about Ro, we expand 
the action Sj(r,r';t) about the points r = R and r' = 
R 0, as 

Si(r,r';t) » S 7 (R ,R ;t)+fiMr-Ro)-«d(r'-Ro), 

(59) 

where a classical particle traveling along the closed path 
7 would leave the point Ro with momentum 7ik T = 
— ^^(Ro, r'; i)/9r'| r / = R and return to Ro after time t 
having momentum ?ik 7 ~ SS^r, R ; t)/dr\ r= n . Sub- 
stitution of Eq. along with Eqs. (|5*H|l and (ff)9"|) into 
Eq. yields 



(</>o|G sc (i)|<M 



x d/2 



2Trih J 



x exp [ / 
where 



&y(Roj Rq; t) 

h 



(60) 



5o(k) 



'ikR,n 



dr|0o(r)|e 



i(k -k)r 



(61) 



In deriving Eq. (|60|l we used -D 7 w D 1 and = v~ 
assuming that the paths, 7 and 7, are close and follow 
the same collision sequence. 

Let us now assume the function |</>o( r )| to be suffi- 
ciently smooth in order for 0o(k), defined by Eq. Q6ip. 
to be sharply peaked about k . Then, the main con- 
tribution to the autocorrelation amplitude in Eq. l(5UI) 
comes from paths 7 with k 7 « k 7 « ko. Thus, in order 
to obtain a leading contribution to (0 O |G sc (i) |0o) a t a 
fixed time t, one can restrict the summation in Eq. 1(60(1 
only to classical periodic orbits 7 passing through a small 
neighborhood of the phase space point (Ro,ko), and 
therefore having the momentum p ss Tiko. The classi- 
cal actions along such periodic orbits can be written as 
5 7 ~ (h 2 kQ/2m)t, and are the same for all 7's. The 
autocorrelation amplitude then reads 



bo\G sc (t)\fo) « e^ fe ' I 1 



2nih 



d/2 



dv\4> Q {v)\) 5> 7 /2 e-^ 



(57) 



(62) 



Equation 1)62(1 is suitable only for predicting the values 
of the autocorrelation function C(t) w |(|</>o|G sc (i)|</>o) | 2 
at times t = i 7 such that there exists at least one classi- 
cal periodic trajectory 7 of period i 7 passing through the 
spatial point Rq with momentum fik^. Due to the narrow 
momentum distribution of the initial wave packet, given 
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by the function <^o(k) in Eq. jSOjl. C{t) decreases by or- 
ders of magnitude as the time t changes to a value such 
that all the classical periodic orbits of period t have mo- 
menta different from hko- Therefore, we come to a con- 
clusion that the time decay of the autocorrelation func- 
tion consists of sequence of sharp peaks centered at the 
times tj. Equation l|62[l is then only suitable for pre- 
dicting the relative strength and time-location of these 
peaks. 

It will be shown below that the probability measure of 
periodic orbits (or D^) decreases exponentially with the 
increase of the number of collisions a classical particle 
undergoes while traveling along these orbits. Therefore, 
the value of the autocorrelation function at a given peak 
at time t is predominantly determined by a subset {7'} of 
the set of all classical periodic trajectories {7} of length 
vt — (hko/m)t, such that the members of the subset have 
the smallest number of scattering events, N(t), possible 
for time t. Indexes Vy = 2N(t) are the same for all 
members of the subset {7'}, yielding 



(63) 



Equation (|63|l allows one to predict the relative mag- 
nitude of peaks of the autocorrelation function C(t) by 
searching (analytically or numerically) for the periodic 
orbits with the smallest number of scattering events dur- 
ing a given time t. We will employ this formula in the 
sequel to calculate the autocorrelation function decay for 
wave packets in various hard-disk and hard-sphere bil- 
liards. 

If properly modified, the above technique is also appli- 
cable for calculation of the peaks of the classical auto- 
correlation function C c \(t), i.e. a fraction of classical tra- 
jectories in a small phase space region around the initial 
location of a classical particle, which return to this re- 
gion after time t. The classical autocorrelation function 
gives the phase space return probability for a classical 
particle described by a phase space distribution function 
rather by the exact coordinates. It characterizes the es- 
cape of classical trajectories from a small neighborhood 
of a chaotic repeller of the system. In mixing chaotic 
systems, such as hard-disk and hard-sphere billiards, the 
autocorrelation function decays exponentially with time, 
C c \(t) ~ e _7clt , with the decay rate 7 c i known as the es- 
cape rate on the repeller 0, E3| • The changes one needs 
to introduce in Eq. (|63[) are apparent. Due to the absence 
of interference in classical mechanics one has to directly 
sum the probabilities of the periodic orbits, Dy, rather 



than the probability amplitudes 



7' j 

/. This leads to 



C cl (t) 



7' 



(64) 



Finally, one needs to specify quantities Dy entering 
the RHS of Eqs. JB3J| and This can be done tak- 

ing into account that the Van Vleck determinant Dy 



equals, up to a normalization factor, the classical proba- 
bility measure of the trajectory 7', e.g. see j2B|. Thus, in 
(i-dimensional space the probability density to find a par- 
ticle at a distance R from the source radiating particles 
with some initial velocity distribution is proportional to 
Furthermore, the probability for a particle to 
undergo a collision with a scatterer in a given direction 
is described by the differential cross section Odiff {&), with 
6 being the scattering angle. Therefore, the probabil- 
ity for the particle to follow a periodic collision sequence 
i]y — {i, j, . . . , q, r, s}, constituting an orbit 7', is 



D 



y 



Cdiff(7T - <t>i) Odiff (tT ~ 4>q) 



R 



d-1 



pd-1 
IXqr 

Odiff (tT - <t>r) Cdiff (vr - 4> s ) 



(65) 



R 



d-l 



R 



where Rij is the center-to-center separation between scat- 
terers i and j, and <pi is the angle of the periodic orbit 
polygon (with the vertices at the scatterer centers) cor- 
responding to the i th vertex. The separation distances 
satisfy an obvious relation, 



R; 



Rn 



Rrs + Rsi — Vt, 



(66) 



which implicitly provides the time dependents to the 
RHS of Eqs. Eg} and 

As we will see below, the expression for the quan- 
tum (classical) autocorrelation function proposed in this 
section, despite its simplicity, accurately predicts times 
and relative magnitudes of the wave packet (distribution 
function) reconstruction peaks. Nevertheless, application 
of more sophisticated techniques, like the one presented 
in the previous chapter, is required if one needs to obtain 
absolute (and not relative) values of the autocorrelation 
function for a wide range of times, including time inter- 
vals between the neighboring peaks. One needs to have 
detailed knowledge of the particle's wave function in or- 
der to predict C(t) for times t other than the peak times, 
i.e. for times that have no phase space period orbits 
of velocity v corresponding to them. This requires the 
construction of the full quantum propagator for a given 
system by methods analogous to the one presented in 
Section II. 



A. Application to studied cases 

We first start with applying the semiclassical technique 
to calculate the peaks of the autocorrelation function for 
the hard-disk scattering systems analyzed in Section II by 
the method of multiple collision expansions, i.e. for the 
two-disk billiard, and three-disk equilateral and isosceles 
billiards. After that we will treat such more complicated 
system as the tree-disk generic billiard, as well as some 
hard-sphere billiards in three dimension. 
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1. Two-disk billiard 

In the case of the two-disk scattering system, see fig El 
there is only one scattering sequence, "1212...", con- 
tributing to a peak of the autocorrelation function, C(t n ), 
at time t n = 2nR/v, with n = 1,2,.... The probability 
weight D ni corresponding to a periodic orbit of length 
vt n , is calculated according to Eq. (|d3|| with d = 2 and 
the semiclassical hard-disk differential cross section 



where 9 = tt for the back scattering. Then, 

2 71 



exp 



(67) 



(68) 



where A^ 2 ^ is the two-disk Lyapunov exponent given by 
Eq. 127|). 

Since the right hand sides of Eqs. lj6l?|) and lj6~4"|) reduce 
to the return probability D n of a single collision sequence, 
we see that the classical autocorrelation function decays 
exactly in the same manner as the quantum one: 



C cl (t n ) ~ exp (-A^t„) . 



(69) 



As we will see later, this similarity is the consequence 
of the fact that the Kolmogorov-Sinai entropy for the 
two-disk periodic orbit is zero [fj, i.e. there is no in- 
formation production in the system since for any time 
t n = 2nR/v there exists only one trajectory leading 
to the wave packet (distribution function) partial recon- 
struction. Thus, the phenomenon of interference between 
different trajectories is absent in the quantum case, re- 
sulting in the same escape rates for classical and semi- 
classical particles. 



2. Three-disk equilateral billiard 

Let us now address the three-disk scattering system, 
with the scatterers centered in the vertices of an equilat- 
eral triangle. Figure 03 shows the autocorrelation peaks 
as a function of time, which were computed numerically 
according to Eqs. (|63[) . I|65|l and l|67|l by summing over 
all collision sequences satisfying Eq. (JHHJ). The dashed 
trend line in the figure represents the exponential de- 
cay with the rate 7^ 3 - ) given by Eq. 1)47(1. The billiard 
is characterized by the disk radius a = 1, and the disk 
center-to-center separation R = 10 4 . The system is iden- 
tical to the billiard considered in Section II, see fig. |3 
and the decay of the autocorrelation peaks is to be com- 
pared with the one presented in fig [3] Note, that the de 
Broglie wavelength A does not enter the sum in Eq. (|63H , 
being contained in a possible prefactor, and is therefore 
not important for determining the relative strength of the 
peaks. Figure |H1 clearly shows that, as in the two-disk 
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FIG. 8: Peaks of the autocorrelation function for the 
equilateral three-disk billiard calculated in accordance with 

Eqs. ((Hill and ffift. The dashed line shows e" 7<3>t de- 

cay, with 7^ 3 - ) given by Eq. (I47II . The radii of the disks consti- 
tuting the billiard equal a = 1, and the disk center-to-center 



separation is R = 

fig- El 



10 . This figure is to be compared with 



case, the semiclassical theory recovers the autocorrela- 
tion function decay rate obtained in Section II by means 
of the multiple collision expansion technique. 

The decay rate 7^ given by Eq. (|47|) for the equi- 
lateral three-disk billiard can be exactly recovered using 
the semiclassical method. In order to calculate 7® one 



needs to sum probability amplitudes 



over all col- 



lision sequences r; 7 ' satisfying Eq. I|66|l . Once again, this 
can be accomplished with the help of the matrix method 
discussed in the previous section. We construct a one- 
collision transition matrix q according to 



1-2 1-3 2 1 2-3 3-1 3-2 
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(70) 



with 



_ / a \V2 



Here, x and 



and 



\fia 
~4R 



1/2 



(71) 



w are the values of the amplitude 
with (Tdiff given by Eq. I(67|l and </> 
The matrix q 



\/crdiff(7r - 4>)/R. 

taking values of and 7r/3 respectively, 
describes a transition due to a single collision event in 
the six-dimensional space spanned by directions (1— »2), 
(l->3), (2->l), (2->3), (3 and (3->2). This matrix 
allows one to express the sum of overlap amplitudes in 
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Eq. I|73|) for times t n = nR/v, with number of collisions 
n = 2, 3, . . ., according to 



(72) 



where the subscript in the RHS denotes that the one- 
one element of the matrix is taken. The autocorrelation 
function at the n th collision is related to the one at the 
(n + l) th collision by 



C(t n + R/v) 
C(t n ) 



(q")i,i 



(73) 



For large number of collisions, n>l, the largest eigen- 
value of the matrix q, equal to x + w, dominates both 
the numerator and the denominator of Eq. (|73[1 , so that 



lim 

t n — > + OC 
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~ R 



C(t n + R/v) 



C{t n ) 

(1/2) 1/2 + (z^/a) 



(x + wf 



1/2' 



exp 



,(3)5 



(74) 



(3) 



where 7 < - 3 * ) is the decay rate given by Eq. (|47|l . 

Finally let us mention that the classical escape rate 7 , 
can be obtained with the help of the transition matrix q, 
given by Eq. I|7U|I . if one redefines x and w by replacing 
the square roots in Eq. I|71[l by the first powers, i.e. 



— and 
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(75) 



Then, 

C cl (t n + R/v) 
C c \(t n ) 
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1 V3 



= exp 
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with the classical escape rate given by 



(3) 
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v , 4i? 
— hi — ==- 
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v , 1.07 R 
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(76) 



(77) 



One can see that the absence of interference in the classi- 
cal case results in faster particle escape from the scatter- 
ing system. The classical escape rate for the three-disk 
equilateral billiard was first obtained by Gaspard and 
Rice 0. 



3. Three-disk isosceles billiard 

In order to complete the comparison of predictions 
of the semiclassical methods with the results of the de- 
tailed binary collision expansion studies, we consider the 
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FIG. 9: Peaks of the autocorrelation function as predicted by 
Eq. 1631 for the isosceles three-disk billiard with a = 5/2. The 

dashed line shows e _7 ° * decay, with 71 3 ' given by Eq. I55L 

while the dotted lines show the trend of the e _A< * decay, 
with A 1 - 2 ' defined by Eq. 1371 . The billiard is parametrized by 
a — 1 and R — 10 4 . This figure is to be compared with fig-Q 



case of the three-disk billiard with the scatterers cen- 
tered in the vertices of an isosceles triangle. Figure 
displays the peaks of the autocorrelation function in 
the system shown in fig. ^ with a = 5/2, a = 1 and 
R = 10 4 . The structure of the decay is twofold. There 
are relatively big recurrences of the wave packet at times 
t n = (2an + l)R/v = (5n + l)R/v with n = 1, 2, 3, . . .. 
The magnitudes of these recurrences follow e~ 7 ° 3> * decay 
(represented by the dashed line) , with the decay rate 7a 
given by Eq. (|55|l . In between any two large neighboring 
peaks the autocorrelation function decays rapidly, ap- 



proximately following e 



-A< 2 >i 



decay (dotted lines), with 



the two-disk Lyapunov exponent defined by Eq. J37J). 
The sequence of autocorrelation function peaks in fig. 
is almost identical to the one in fig. 

We will now use the matrix method to derive the au- 

(3) 

tocorrelation function envelope decay rate 7 Q directly 
from Eq. H63|) . As it was mentioned above the sum in 
the RHS of Eq. (|rJ3*|> goes only over such periodic or- 
bits 7' that comprise the smallest number of scatter- 
ing events possible for a periodic trajectory of length 
vt. This is because the overlap amplitudes y ' D^* cor- 
responding to trajectories with the longest mean free 
path are given by products of the smallest number of 
\f^dis/R *C 1 terms, and therefore are the dominant 
ones. In the case of the three-disk isosceles billiard with 
a > 1 these long mean free path trajectories are the ones 
that pass through the scatter "3" every second collision, 
see fig. 0] This amounts to the collision sequence "132" 
at time t\ = (2a + l)R/v, collision sequences "13132" 
and "13232" at time t 2 = (4a + l)R/v, etc. In gen- 
eral there are 2 n_1 different periodic orbits of length 
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vt n = (2na + l)R contributing to the autocorrelation 
function peak C(t n ). 

The relative probability weights of the above trajecto- 
ries are determined according to Eq. (|65|) with the dif- 
ferential cross section given by Eq. JB7J. Thus, collisions 
(1 — ► 3 — ► 1) and (2-^3^2) are described by the cross 
section (Jdiff3 = a /2, while (1^3—^2) and (2— >3— > 1) by 
ffdiffa = (a/2)cos(0 3 /2) = (a/2)y/l - l/4a 2 , where </> 3 
denotes the angle of the triangle corresponding to vertex 
"3" , see fig. 0J The first and the last collisions of every 
periodic trajectory deflects the moving particle by the 
angle 6 = it — <pi = tv/2 + <j>$/2, and is described by the 
cross section Cdiffi = 0"diff2 = (a/2 3 / 2 )-^/l + l/2a. Here 
<pi and </>2(=</>i) are the other two angles of the isosceles 
triangle satisfying 2(pi + tfi 3 — n. Thus, the sum of over- 
lap amplitudes at time t n — (2na + l)R/v, corresponding 
to periodic orbits of 2n + 1 collisions, can be written as 
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where 
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with 
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and 
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(80) 



Repeating the arguments used in the case of the three- 
disk equilateral billiard we find that 



C(t r 



x 2 (x + w) 




(81) 



since \J x(x + w) is the largest eigenvalue of the matrix q. 
The time interval between any two successive large peaks 
of the autocorrelation function is t n +i — t n = 2aR/v, so 
that for n ^ 1 we get 



C{t n+1 ) w C{t n ) exp - 7 ^)(t 



with the decay rate given by Eq. I|55|l . Once again we 
observe strong agreement between the prediction of the 
semiclassical analysis and the results of the diffraction 
regime approximation obtained in Section II. 

(3) 

The classical decay rate 7^ a for the three-disk isosce- 
les billiard can be obtained by the following modification 
of the elements of the transition matrix q: 
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Then, 
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(84) 



with the classical decay rate 
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(85) 



Comparison of Eqs. (|55fl and (|85fl shows that, as in the 
case of the three-disk equilateral billiard, classical parti- 
cle escape in the isosceles billiard takes place at higher 
rate than the corresponding quantum process. 



B. More billiards 

As we have seen, the semiclassical method presented 
here is suitable for predicting the main features of the au- 
tocorrelation decay in hard-disk scattering systems. We 
now apply the method to problems which could not be 
easily treated by the technique of explicit calculation of 
scattering resonances used in the Section II. 



1. Generic three-disk billiard 

The first system we address is a three-disk billiard of 
the most general type: the disks of radii a are centered 
in the vertices of a triangle of unequal sides R, aR and 
f3R, where for concreteness we take [3 > 1, a. In order 
to visualize the system one should consider the three- 
disk isosceles billiard shown in fig. and elongate the 
side "13" of the triangle from its original length aR to 
the new length f3R. It is a formidable problem to cal- 
culate the scattering resonances (and the corresponding 
residues) governing the time evolution of a wave packet 
for such a system. This makes the multiple collision ex- 
pansion technique, used to study the three-disk equilat- 
eral and isosceles billiards, inefficient in the case of the 
generic three-disk scattering system. On the other hand 
the semiclassical approach of this section is quite easy to 
implement. It allows one to extract such important infor- 
mation about the autocorrelation function as the relative 
strength of the revival peaks and the overall envelope de- 
cay rate. 

The trajectories with the longest free flight path are 
the ones that bounce most of their time between disks 
"1" and "3" with the largest center-to-center separation 
PR. In the limit of a large number of collisions, n ^> 1, 
there is a single periodic collision sequence "1313. . . 132" 
resulting in a relatively strong wave packet recurrence at 
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times t n = [1 + a + (n — 2)(3}R/v. In accordance with 
Eqs. !p5)l. and $67$ we have 



C{t n ) 



( a \ ™ 2 a { <p 2 \ a (fa 

\m) w cos {t) 2^r cos {t 



exp 



(86) 



where 02 and 4>3 are the triangle angles at vertices "2" 
and "3" respectively, and 



(87) 



is the two-disk Lyapunov exponent, see Eq. (|37|l . corre- 
sponding to disks "1" and "3" . 



of neighboring peaks, so that only the overall envelope 

. (2) 

decay C(t) ~ e p can be resolved. Thus, the simple 
semiclassical approach of the section allows one to pre- 
dict main features of the time-dependent autocorrelation 
function for wave packets in arbitrary shaped three-disk 
billiards. 

Finally, let us compare the decays of the wave packet 
autocorrelation function in three-disk scattering systems 
of three possible types: (i) equilateral (a — (3 = 1), (ii) 
isosceles (a = |3 / 1) and (iii) generic (a / |3 / 1) 
three-disk billiards. Let also the triangles, constituted by 
the disks centers, have approximately equal side lengths, 
and differ only by the number of symmetries. Thus, in 
the cases (ii) and (iii) both a and (3 are close (but not 
identical) to unity. Table^represents the autocorrelation 




FIG. 10: Peaks of the autocorrelation function for the three- 
disk billiard with a = %/2, (3 = 2, a = 1 and R = 10 4 . Solid 

line represents e X < 3 * decay, with \p calculated according 

to Eq. 1871 : dashed line corresponds to e~ A< * decay, with 
given by Eq. lEiTl . 



Figure EH shows peaks of the wave packet autocorre- 
lation function for the three-disk billiard with a = 
(3 = 2, a = 1 and R = 10 4 . The peaks fall inside a narrow 
cone. Magnitudes of the relatively strongest recurrence 



peaks decay exponentially with time as e 



with the 



decay rate A^ predicted by Eq. (J57J). The trend of this 
exponential decay is shown by the solid line. The dashed 

line represents the exponential decay e _A< '* due to the 
shortest two-disk periodic orbit in the system. The value 
of A^ 2 ) is calculated in accordance with Eq. 1)3 7|) . and 
represents the fastest decay rate in the billiard. 

As shown in Section II, the peaks of the autocorrela- 
tion function have significant width which increases with 
time, see figs. |3 |S] and \7\ Therefore, if the autocorre- 
lation function peaks are closely spaced, as in fig. 1101 
the peak broadening ultimately results in overlapping 



Disk billiard 


Decay exponent 


Equilateral 

a = (3 = l 


v , 0.54 R 
— m 

R a 


Isosceles 

a = (3^1 


v , 1.04 i? 
— In 

R a 


Generic 
a^/3^1 


v , 2R 
— In — 
R a 



TABLE I: Autocorrelation function decay exponents for 
three-disk billiards of different symmetries. All billiard are 
based on almost equilateral triangles, i.e. a and f3 in the 
second and third row of the table are close to unity. 

function decay exponents for the above-mentioned three- 
disk billiards calculated in accordance with Eqs. I|47|l . 
(|55|l and l|87|) . One can see that the decay rate increases 
by approximately (v / R) ln2 as the number of equal sides 
in the three-disk billiard decreases by one. This difference 
in the decay rates is approximately equal to twice the 
difference in KS-entropies per unit time of corresponding 
billiards, and will be discussed in more details in the 
sequel. 

Presence of symmetries in a scattering system increases 
the number of periodic trajectories of a given length, and 
thus enhances interference effects. It is due to the inter- 
ference that strong wave packet reconstruction peaks oc- 
cur resulting in slower envelope decay of the wave packet 
autocorrelation function. 



2. Hard-sphere billiards in three dimensions 

The above semiclassical approach can also be applied 
to calculate the autocorrelation function decay for wave 
packets in three-dimensional hard-sphere billiards. Here 
we briefly derive the autocorrelation function decay rates 
for the following three systems: (i) the two-sphere, (ii) 
the three-sphere equilateral, and (iii) the four-sphere 
tetrahedral Eal billiards. 
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The wave packet (phase space distribution) partial re- 
construction peaks of the semiclassical (classical) auto- 
correlation function are given by Eq. (|63(1 (Eq. I|fl4)l). with 
the collision sequences 7j 7 ' satisfying Eq. (|66|) . The prob- 
ability weight Dy of the periodic orbits 7' are now deter- 
mined from Eq. Ij65(l with d = 3, and with the differential 
cross section 



of the three-dimensional transition matrix 



Cdiff(60 



4 ' 



where a stands for the radius of the hard-sphere scatterer. 
It will become clear below that the independence of the 
differential cross section ddiff on the scattering angle 8 
significantly simplifies the calculation of the decay rates 
of C(t) and C c i(i). 



Let us start with analyzing the autocorrelation func- 
tion peaks for a wave packet moving in a two-sphere bil- 
liard consisting of two hard spheres, "1" and "2" , of ra- 
dius a, separated by a distance R^> a, e.g. see fig. |21 A 
wave packet is initially located on the line connecting the 
sphere centers and moves toward one of the spheres. As 
in the two-disk billiard case, there exists only one periodic 
collision sequence, "1212. . . " , contributing to the partial 
wave packet (phase space distribution) reconstruction at 
times t n — 2nR/v, with n = 1,2,.... Hence, 



C(t n ) ~ C c \(t n ) 



\4R 2 J 



2n 



exp 



where 



7 (2) -7 (2) -2\W-— In — 

'3D Tcl,3D ZA R ■ 



(2) f 
"73D l n 



2R 



(89) 



(90) 



is the sum of the two positive Lyapunov exponents, 
A^ 2 ' = (v/R)ln(2R/a), of the two-sphere periodic orbit. 
As in the two-disk billiard case, the absence of the inter- 
ference of different trajectories results in equality of the 
semiclassical and classical decay rates. 



The equilateral three-sphere billiard in three dimen- 
sions consists of three spheres of radius a placed in the 
vertices of an equilateral triangle with side R >■ a, 
e.g. see fig. 0J In this system, the peaks of the au- 
tocorrelation function occur at times t n — nR/v, with 
n = 2,3,4,... counting the number of collisions. The 
probability amplitude of a closed paths 7', given by 
y/Dy(t n ) = (a/2R) n , is only a function of time t n and 
does not depend on the details of the particular collision 
sequence r) 7 > = . . . , q, r, s} constituting the orbit 7'. 
The number M(t n ) of all collision sequences r]j> satis- 
fying Eq. l|66f) for t = t n , and therefore the number of 
interfering trajectories, can be calculated with the help 
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(91) 



in accordance with 



M(t n ) = (q£ D ) M . 



(92) 



Here, as above, the subscript "1,1" denotes that the one- 
one element of the matrix is taken. Then, the magnitude 
of the autocorrelation function peak at time t n is given 
by C(t n ) - [M{t n )(a/2R) n ] 2 . For long times, n > 1, the 
decrease of the peak strength due to one collision can be 
written as 



C{t n + R/v) 
C(t n ) 



(qSD)x,i 



(—)' 

\2RJ 



\R 



(93) 

Here we used that for n 3> 1 the ratio of the matrix ele- 
ments in the last equation is given by the largest eigen- 
value of the matrix q3D, which is equal to 2. Therefore, 
we have 



C(t n ) ~ exp (-730 t^j , 



where 



73D R a 



2v , R 



(94) 



(95) 



is the autocorrelation function decay rate for the three- 
sphere equilateral billiard in three dimensions. 

The classical decay rate is calculated in the analo- 
gous way. According to Eq. (|64|l we write C c \(t n ) ~ 
M(t n ){a/2R) 2n . Then, 



n+l\ 



c cl (t n + R/ v ) (q' 3 D 



Cci(t n ) (qJ D ) M \2RJ 

and consequently 

C c \(t n ) - exp f-7cf,3D *n) 



V2R 



where 



(3) _2v V2R 

7cl,3D ~ R ln a 



(96) 



(97) 



(98) 



is the classical autocorrelation function decay rate of the 
three-sphere equilateral billiard. Below we will see that 

(3) (3) 

the difference between 7^ 3D and 7 3D is equal to the topo- 
logical entropy per unit time (which in this particular 
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FIG. 11: The tetrahedral four-sphere billiard. The wave 
packet is initially placed between spheres "1" and "2", with 
its average momentum directed toward sphere "2" . 



system coincides with the KS-entropy per unit time) of 
the chaotic repeller of the classical system. 

Finally, we consider a substantially three-dimensional 
scattering system - a tetrahedral billiard - where four 
spheres of radius a are placed in the vertices of a pyramid 
build of four equilateral triangles of sides R a. The 
wave packet starts on a line connecting two disks labeled 
by "1" and "2" , see fig. El 

The calculation of the autocorrelation function decay 
exponent for this billiard proceeds in close analogy with 
the three-sphere case considered above. The transition 
matrix qjD is now a 12 x 12 matrix connecting directions 
(1 - 2), (1 -> 3), (1 -> 4), (2 - 1), (2 - 3), (2 -> 4), 
(3 -» 1), (3 -» 2), (3 -> 4), (4 -> 1), (4 -> 2) and 
(4 — > 3). The structure of the matrix is similar to the 
one given by Eq. I|91[) for the case of the three-sphere 
billiard; the largest eigenvalue is now equal to 3. As 
before the number of periodic trajectories contributing to 
the autocorrelation peak at time t n = nR/v, n = 2, 3, . . ., 
is given by Eq. H92|l with q3D being the 12 x 12 transition 
matrix of the four-sphere tetrahedral billiard. Then, the 
decrease of the autocorrelation function peak strength 
due to one collision is given by 



C(t n + R/v) fi^) 1,1 



C{t n ) 



leading to 



(—)' 

\2RJ 



3n 
2R 



(99) 



C(t n ) ~ exp f-7gn tn) 



with the four-sphere escape decay rate given by 



(4) 

7 3D 



2v , 2R 
— In . 

R 3a 



(100) 



(101) 



The classical autocorrelation function decay rate for 
the four-sphere tetrahedral billiard is obtained straight- 
forwardly. According to Eq. (|64|l we write 



c cl (t n + R/ v ) (q^M 



C c \(t n ) 
This results to 

with 



(q« D ) M \2R) 

C c \{t n ) ~ exp (-7cL3D O 




7 C 1,3D 



2« , 2i? 
— In . 

-R v3a 

(4) 



(104) 



The semiclassical decay rate, J^, is again found to be 

slower than the classical one, 7cf 3D - As we will show be- 
low the difference between the two decay rates is given by 
the topological entropy per unit time of the four-sphere 
tetrahedral billiard repeller. 
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Four sphere tetrahedral billiard 
Three sphere equilateral billiard 
Two sphere billiard 




FIG. 12: Peaks of the wave packet autocorrelation function 
for the three scattering systems: two-sphere (circles), three- 
sphere equilateral (triangles) and four-sphere tetrahedral (di- 
amonds) billiards. In all cases the radius of the sphere scat- 
terers is a = 1, and the sphere center-to-center separation 
is R = 10 4 . The semiclassical decay rates for these billiards 
are calculated according to Eqs. I|9()|l . I|95^ and IjlOljl . and 
presented in the figure by the dotted, dashed and solid lines 
respectively. 

Figure El shows relative magnitude of peaks of the 
wave packet autocorrelation function for the three scat- 
tering systems considered above: two-sphere, three- 
sphere equilateral and four-sphere tetrahedral billiards. 
The radius of the sphere scatterers is a = 1, and the 
sphere center-to-center separation is R — 10 4 . The semi- 
classical autocorrelation function decay rates for these 
billiards are calculated according to Eqs. (|90|l . (|95|l and 
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pOip . and presented in the figure by the dotted, dashed 
and solid lines respectively. The figure illustrates the re- 
duction of escape of a semiclassical particle from a billiard 
as additional scatterers are added to the latter. Below 
we clarify the connection of the rate of this escape to 
classical properties of the system's chaotic repeller. 

C. Classical and semiclassical escape rates 

We would now like to comment on the connection of 
the escapes rates in classical and semiclassical open bil- 
liards to such properties of classical chaotic systems as 
the mean Lyapunov exponents, the Kolmogorov-Sinai 
(KS) and topological entropy. 

The classical escape rate 7 c i, which is equivalent to 
the decay rate of the classical autocorrelation function, 
is known to equal a difference of the stretching and ran- 
domization rates in a chaotic system |30|, |3jJ, |38| , 

7ci= ( 105 ) 

A;>0 

Here, the sum, going over all mean positive Lyapunov 
exponents A<j , represents the rate of the local exponential 
stretching of an initial particle phase space distribution, 
and /iks is the KS-entropy per unit time of the system 
characterizing the rate at which the phase space distri- 
bution gets randomized over the chaotic repeller of the 
scattering system. It was first pointed out by Gaspard 
and Rice Q that the semiclassical escape rate 7 for the 
three hard-disk scattering system in two dimensions is 
given by 

7«A-2/i KS , (106) 

where A is the mean positive Lyapunov exponent of the 
corresponding chaotic repeller. Comparison of Eqs. 1)105(1 
and (|106Jl shows that the semiclassical escape rate is 
never greater than the escape rate in the counterpart 
classical system, which is consistent with our observa- 
tions earlier in this section. 

It is interesting to note that the equality in Eq. I|10tj|) 
can be extended (and made exact!) for the three- 
dimensional hard-sphere billiards considered in this sec- 
tion, i.e. for the two-sphere, three-sphere equilateral 
and four-sphere tetrahedral scattering systems. More- 
over, Eq. (|106f) can be easily derived for these billiards, 
since the KS-entropy of the repellers in these systems 
has an especially simple form. Indeed, the independence 
of the hard-sphere differential cross section on the scat- 
tering angle, evident from Eq. (|88(l . and the equality 
of lengths of all free flight segments composing the sys- 
tem's periodic orbit, result in equivalence of probability 
weights for all trajectories of given length constituting 
the chaotic repeller, see Eq. I|65|l . Since the probability 
measure is the same for all periodic trajectories involving 
a given number of collisions, the construction of the KS- 
entropy is identical to the construction of the topological 



entropy, e.g. see j^- Therefore, for the particular three- 
dimensional billiards of this section, the KS-entropy per 
unit time, /iks, of the repeller is simply equal to the 
topological entropy per unit time, h top - The later is the 
rate of the exponential growth with time t of the number 
of different collision sequences, M(t), a classical particle 
moving on the repeller can possibly undergo: 

Af(t)~exp(fct„pt). (107) 

For the two-sphere billiard there is only one periodic 
collision sequence a particle staying on the repeller can 
follow. Thus, M (t) = 1 and 

hill = °- ( 108 ) 

In the three-sphere equilateral billiard the number of pos- 
sible trajectories multiplies by two at every collision, so 
that M(t) = 2 vt / R , leading to 

hfX = ^\n2. (109) 

In the same way we conclude that in the four-sphere 
tetrahedral billiard the number of possible collision se- 
quences is M(t) = 3 vt ' R , so that the topological entropy 
per unit time reads 

^ = 1^3. (110) 

Every trajectory 7' of the three-dimensional hard- 
sphere billiard repeller has two positive Lyapunov ex- 
ponents, Ai(7') and A2(7'), corresponding to perturba- 
tions of the particle's initial conditions in two direction 
perpendicular to the trajectory. It is straightforward to 
show that the sum of the two Lyapunov exponent, for a 
dilute hard-sphere scattering system, R ^> a, does not 
depend on the details of the trajectory 7', and is given 
by [H 

2v 2R , . 

Ai + A 2 = — In— . Ill 
R a 

The last equality is obvious for the two-sphere peri- 
odic orbit, since in that case Ai = A2 = A^ 2 ^ = 
(v/R)\n(2R/a). In the case of a general hard-sphere 
dilute scattering system the method of curvature radii 
[il^. l30| can be used to prove Eq. (|111|) . 

Looking at the expressions for the classical autocorre- 
lation function decay rates in the three-dimensional bil- 
liards studied above, Eqs. l|9"8|) and (|104fl . we notice 
that the following relation holds 

tSd = Ai + A 2 - h®, (112) 

with j — 2, 3, 4. One can also verify that the semiclassical 
decay rates, given by Eqs. lj9"U|) . l|9"5|l and (|101fl . satisfy 

7$ = Ai + A 2 - 2h<£ p . (113) 
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Equations i|112fl and i|113|) are the three-dimensional ver- 
sion of Eqs. I|1U5|) and (|106f) respectively. 

The appearance of the factor of 2 in the expression for 
the semiclassical decay rate is now apparent: the strength 
of the wave packet partial reconstruction peaks, C(t), is 
proportional to M 2 (t) due to interference of different pe- 
riodic orbits, while in the classical case C c \(t) is given 
merely by the sum of orbit probabilities, and is propor- 
tional to M(t). One can also see now that the difference 
of the two escape rates is given by the topological entropy 
(which, for the three-dimensional billiards in question, is 
equal to the KS-entropy) of the chaotic repeller of the 
classical system, and is therefore determined by underly- 
ing classical dynamics. 




r 



Ay(t) 



FIG. 13: Schematic picture of a classical periodic orbit 7' to- 
gether with an exponentially spreading tube of infinitesimally 
close trajectories. The initial cross section of the tube, j4 7 < (0), 
gets magnified to the value j4y(t), given by Eq. 11141 . after 
time t equal to the period of the orbit 7'. 



Finally, we would like to further clarify the appearance 
of the factor of 2 in the expression for the semiclassical 
escape rate, Eq. ()llr>(l . by means of the following sim- 
ple argument. Consider a classical periodic orbit 7' of 
the period t, schematically depicted in fig. ^2 passing 
through the phase space point around which the initial 
wave packet is localized. In a general chaotic system 
the number of such periodic orbits M(t) grows exponen- 
tially with the period t, and the rate of this growth is 
given by the topological entropy per unit time, /itop, see 
Eq. (|1(J7[I . In order to determine the probability weight 
Dji(t) of the orbit 7', one needs to consider a bundle of 
trajectories which stay infinitesimally close to 7'. These 
trajectories form an exponentially thickening tube, see 
fig- El with the initial and final cross sections, ^4 7 '(0) 
and Aji (t) respectively, related by 



A Y (t)= Ay (0)exp\ W)* 



(114) 



,A J ( 7 ')>0 



Here the sum goes over all positive Lyapunov exponents 
Xj(j') of the periodic orbit 7'. Then, the probability for 
a classical particle taken at random from the trajectory 
bundle to stay on the periodic orbit 7' is proportional to 



the ratio Ay/(0)/Ay/(i), so that 



LV(*)~exp(- Y W)* 

M7')>0 



(115) 



Equation l|115|) together with (|107fl leads to the expres- 
sion for the semiclassical escape rate in the case when the 
sum over positive Lyapunov exponents, X^Aj'CtOi does 
not significantly depend on the details of the orbit 7', 
and can be replaced by its average value ^ <\? . Then, 
the sum over orbits {7'} in Eq. (|rJ5|> is simply equal to 
the product of the number of these orbits, M{t), given by 
Eq. (|107fl . and the probability amplitude, -J D(t), iden- 
tical for all the orbits and calculated in accordance with 
Eq. i|115l) . The substitution yields 



C{t) ~ (M(t)V^) 



exp 




2h 



top 



. (U6) 

Thus, we recover the expression for the semiclassical au- 
tocorrelation function decay rate, Eq. (|113J) . obtained for 
a number of hard-sphere billiards in three spatial dimen- 
sions. 

The purpose of the above oversimplified derivation is 
only to illustrate the origin of the intimate relation be- 
tween the quantum escape rate and such properties of 
the underlying classical chaotic system as the Lyapunov 
exponents and topological (and/or KS) entropy. The de- 
tails of this relation as well as the limits of its applicabil- 
ity are yet to be investigated. 



IV. SUMMARY AND CONCLUSIONS 

In the first part of this paper we used the technique of 
multiple collision expansions to construct the quantum 
propagator for a particle with the de Broglie wavelength, 
A, traveling in an array of hard-disk scatterers of radius 
a>A. The scattering system was assumed to be so di- 
lute that the typical scatterer separation R satisfied the 
condition R^> a 2 / A. The quantum propagator was used 
to analytically calculate the time-dependent autocorre- 
lation function for a wave packet, initially localized in 
both position and momentum spaces, evolving in open 
two- and three-disk billiard systems. It was found that 
the autocorrelation function exhibits a sequence of sharp 
peaks at times multiple to periods of classical phase space 
periodic orbits of the billiards. These peaks correspond 
to partial reconstructions of the initial wave packet in 
the course of its time evolution. The envelope of the 
correlation function decays exponentially with time after 
one or two particle-disk collisions; the exponential decay 
lasts for some a/X » 1 scattering events. Our calcula- 
tions recovered the autocorrelation function decay rate, 
first obtained in reference Q , for a particle moving in the 
three-disk equilateral billiard, and predicted the detailed 
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structure and the decay rate of the autocorrelation func- 
tion for the more complicated three-disk isosceles billiard. 

In the second part of this paper the method of the 
semiclassical Van Vleck propagator was utilized to de- 
rive a simple expression for the relative magnitude of 
the wave packet partial reconstruction peaks. Although 
fails to describe the full time dependence of the auto- 
correlation function, this method allows one to calculate, 
analytically or numerically, the decay rates of the au- 
tocorrelation function envelope (also known as particle 
escape rates) for much more complicated hard-disk and 
hard-sphere scattering systems. Here we used it to ana- 
lyze the three-disk generic billiard, the two-sphere, three- 
sphere equilateral and three-sphere tetrahedral scattering 
systems. 

A straightforward modification of the semiclassical 
method allowed us to construct the expression for the 
peaks of the classical autocorrelation function, and to 
compare classical and semiclassical particle escape rates 
in various open hard-disk and hard-sphere billiards. The 
semiclassical escape rate in three-dimensional scattering 
systems was shown to be given by the difference of the 
sum of the two positive Lyapunov exponents and twice 
the topological entropy (in this case equal to the KS- 
entropy) per unit time of the underlying classical re- 
peller. Thus, the semiclassical escape rate is never grater 
then the classical one, and the difference between the two 
equals to the topological (or KS) entropy per unit time 
of the classical system. This result is consistent with the 
earlier findings of reference for the case of the two- 
dimensional hard-disk billiards, and therefore strength- 
ens the connection between classical and quantum chaos. 

An interesting seeming "paradox" arises if one com- 
pares the expressions for the classical and semiclassical 
escape rates, e.g. see Eqs. I|112[) and (|113|) . Indeed, ac- 
cording to the Correspondence Principle it is expected 
that in the limit of the de Broglie wavelength going to 
zero the semiclassical escape rate should match its clas- 
sical counterpart. This seems to be incompatible with 
the results of this paper: the semiclassical escape rate is 
expressed in terms of classical quantities only (and does 
not depend on the de Broglie wavelength), and (in the 
case of hard-sphere billiards) is smaller than the classical 
one by the value of the topological entropy of the system. 
In fact, the difference of the two escape rates originates 
from noncommutivity of the order in which the infinite 
time limit and the classical (de Broglie wavelength going 
to zero) limit are taken. This question has been carefully 



studied by Barra and Gaspard for the case of quantum 
graphs [33, as well as discussed in Ref. 40]. 
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APPENDIX A: CIRCULAR WAVE PACKET 
EXPANSION 

Here we derive the expansion of the circular wave 
packet given by Eqs. (|17|l and l|18|l . The function 
Xi(k, fcn), defined in Eq. Q18[l. can be put in the integral 
form 1271 



Xi(k,k ) = / drrJi(kr)Ji(k r) 



(Al) 



Then, 



/— — OO 



/ rfrr V J ; (fcr)J ; (fc r)e li 

^ r*TJ (|k-k |r), 
a Jo 



(A2) 



where the "summation theorem" for zeroth order Bessel 
function [^3 was used. Doing the simple integral we end 
up with the equality 



£ xHMoK^o) = 2^ J ° ( ' k \ k f ) . (A3) 



[1] F. Haake, Quantum singnatures of chaos, 2nd Ed. 

(Springer, Berlin, 2000). 
[2] H. J. Stockmann, Quantum Chaos: An Introduction, 

(Cambridge U. Press, Cambridge, 1999). 
[3] M. V. Berry, Ann. Phys. 131, 163 (1980). 
[4] M. Sieber and K. Richter, Phys. Script. T90, 128 (2001). 
[5] H. Primack and U. Smilansky, Phys. Rev. Lett. 74, 4831 

(1995); Phys. Rep. 327, 1 (2000). 



[6] P. Gaspard and S. A. Rice, J. Chem Phys. 90, 2225 
(1989). 

[7] P. Gaspard and S. A. Rice, J. Chem Phys. 90, 2242 
(1989). 

[8] P. Gaspard and S. A. Rice, J. Chem Phys. 90, 2255 
(1989). 

[9] P. Gaspard, D. Alonso, T. Okuda, and K. Nakamura, 
Phys. Rev. E 50, 2591 (1994). 



24 



[10] A. Wirzba, Chaos 2, 77 (1992); A. Wirzba and M. 

Henseler, J. Phys. A 31, 2155 (1998); A. Wirzba, Phys. 

Rep. 309, 1 (1999). 
[11] A. Wirzba, "The 1-disk propagator in scalar problems 

and elastodynamics" (2001), unpublished. 
[12] P. Cvitanovic and B. Eckhardt, Phys. Rev. Lett. 63, 823 

(1989). 

[13] E. Doran and U. Smilansky, Nonlinearity 5, 1055 (1992); 

R. Bliimel, B. Dietz, C. Jung, and U. Smilansky, J. Phys. 

A 25, 1483 (1992). 
[14] A. Peres, Ann. Phys. 129, 33 (1980). 
[15] E. J. Heller, E. B. Stechel and M. J. Davis, J. Chem. 

Phys. 73, 4720 (1980). 
[16] E. J. Heller, Les Houches, Session LII, 1989, Chaos and 

Quantum Physics, M.-J. Giannoni, A. Voros, and J. Zinn- 

Justin, eds. (North Holland, Amsterdam, 1991). 
[17] S. Tomsovic and J. H. Lefebvre, Phys. Rev. Lett. 79, 

3629, (1997). 

[18] R. W. Robinett, Phys. Rep. 392, 1-119 (2004). 
[19] A. Goussev and J. R. Dorfman, Phys. Rev. E 71, 026225 
(2005). 

[20] A. Peres, Phys. Rev. A 30, 1610, (1984). 

[21] T. Prosen, T. H. Seligman a nd M. Znidaric, "Theory of 
quantum Loschmidt echoes", quant-ph/0304104 

[22] R. A. Jalabert and H. M. Pastawski, Phys. Rev. Lett., 
86, 2490, (2001); F. M. Cucchietti, H. M. Pastawski, 
and D. A. Wisniacki, Phys. Rev. E65 045206, (2002); F. 
M. Cucchietti, C. H. Lewenkopf, E. R. Mucciolo, H. M. 
Pastawski, and R. O. Vallejos, Phys. Rev. E 65, 046209 
(2002); F. M. Cucchietti, H. M. Pastawski, and R. A. 
Jalabert, Phys. Rev. B 70, 035311 (2004). 

[23] K. M. Watson, Phys. Rev. 105, 1388 (1957). 

[24] C. J. Joachain, Quantum Collision Theory, 3rd Ed. 
(North Holland, Amsterdam, 1983). 

[25] S. Correia, J. Phys. A 34, 9349 (2001); 
S. Correia, PhD Dissertation, Universite 
Louis Pasteur, Strasbourg (2000), webpage 
|http://lptl.u -strasbg.fr/correia/main_these.ps.gz] 

[26] G. B. Arfken and H. J. Weber, Mathematical methods for 
physicists, 4th Ed., (Academic Press, 1995), pp. 671-674. 

[27] I. S. Gradshteyn and I. M. Ryzhik, Table of integrals, 
series, and products, 6th Ed., (Academic Press, 2000). 

[28] H. M. Nussenzveig and W. J. Wiscombe, Phys. Rev. A 
43, 2093 (1991). 

[29] A. Goussev, PhD Dissertation, University of Maryland, 
College Park (2005). 

[30] P. Gaspard, Chaos, scattering theory and statistical me- 
chanics, (Cambridge University Press, Cambridge, 1998). 



[31] E. Bogomolny, Nonlinearity 5, 805 (1992); Chaos 2, 5 

(1992) . 

[32] O. Agam and S. Fishman, J. Phys. A 29, 2013 (1996). 

[33] P. Cvitanovic, Phys. Rev. Lett. 61, 2729 (1988). 

[34] M. C. Gutzwiller, Chaos in Classical and Quantum Me- 
chanics, (Springer- Verlag, New York, 1990); M. Brack 
and R. K. Bhaduri, Semiclassical Physics, (Addison- 
Wesley, Reading, 1997). 

[35] S. Tomsovic and E. J. Heller, Phys. Rev. E 47, 282 

(1993) . 

[36] A. Peres, Quantum Theory: Concepts and Methods, 

(Kluwer Acad. Publ., Dordrecht, 1995). 
[37] H. Kantz and P. Grassberger, Physica D 17, 75 (1985); 

J.-P. Eckmann and D. Ruelle, Rev. Mod. Phys. 57, 617 

(1985). 

[38] J. R. Dorfman, An Introduction to Chaos in Nonequilib- 
rium Statistical Mechanics, (Cambridge University Press, 
Cambridge, 1999). 

[39] F. Barra and P. Gaspard, Phys. Rev. E 65, 016205 
(2001). 

[40] P. Gaspard, NATO ASI Series, Dynamics: Models and 
Kinetic Methods for Nonequilibrium Many-Body Sys- 
tems, J. Karkheck, ed. (Kluwer Academic Publishers, 
Dordrecht, 2000) pp. 425-456. 

[41] We have chosen an initial wave packet that allows us to 
close the integration contour. This would not be possible 
for a Gaussian initial wave packet, and one would need 
to employ different methods to calculate the integral in 
Eq. 

[42] The decay of the free streaming part of the autocorrela- 
tion overlap, Qo(i), results from the spatial separation of 
the wave packets at times and t, and is predominantly 
Gaussian. Reference [2^1 estimates the magnitude of this 
overlap to be thousands orders of magnitude smaller than 
the one of the scattering part of the autocorrelation over- 
lap, fis(i), for the set of parameters used to illustrate 
the C(t) decay throughout this paper, i.e. a = a = 1, 
R= 10 4 and %= 10 -2 . 

[43] The situation is similar to the one observed in the short 
time spreading of small Gaussian wave packets in the 
Lorentz gas systems, e.g. see [Tfll ]. 

[44] The small parameter in the above expansion is the ratio 
of spatial distances of r, r' from Ro to the total distance 
travelled by a classical particle in the periodic orbit. 

[45] The tetrahedron is a triangular pyramid having congru- 
ent equilateral triangles for each of its faces. 



